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We present a rigorous derivation of classical molecular dynamics (MD) from quan- 
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1 Introduction 

A basic mathematical formulation of the passage from quantum to classical mechanics, following 
the ideas of Eugene Wigner |Wig32| , is due to Lions and Paul [LP93j and Gerard [Ger91aj : if 



{^ e } is a sequence of solutions to a semiclassicaliy scaled Schrodinger equation with smooth 
potential, then the sequence of corresponding Wigner transforms converges (up to subsequences) 
to a solution of the Liouville equation, i.e. the transport equation of the underlying classical 
dynamics. For a closely related mathematical approach to semiclassical limits, which goes back 
to Egorov and is based on Weyl quantization and Moyal calculus, see e.g. |Rob87l IMar02] . 

Due to the reliance on smooth potentials, and in particular on the existence and uniqueness 
of trajectories of the classical dynamics, these results are not directly applicable when one 
tries to derive classical molecular dynamics (MD) from Born-Oppenheimer quantum molecular 
dynamics (QMD). By the latter, one means quantum dynamics of the molecule's atomic nuclei 
in the exact non-relativistic Born-Oppenheimer potential energy surface given by the ground 
state eigenvalue of the electronic Hamiltonian with Coulomb interactions. The limit where 
the natural small parameter in QMD, the ratio of electronic to nuclear mass m e /m n =: e 2 , 
tends to zero, has the structure of a semiclassical limit. (The physical value of this parameter 
is ~ 1/2000 for hydrogen, and even less for the other atoms.) However, the potential energy 
surface of QMD is not even continuous, because it always contains Coulomb singularities due to 
nuclei-nuclei repulsion; in addition it can have cone-type singularities at electronic eigenvalue 
crossings. 

Here we present a rigorous derivation of MD from QMD in the limit of small mass ratio 
that is applicable to the exact Born-Oppenheimer potential energy surface with Coulomb in- 
teractions. Our result is valid away from eigenvalue crossings. This is done by extending the 
approach of Lions and Paul |LP93j to an appropriate class of non-smooth potentials. The 
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main technical novelty is a non-concentration estimate on the set of Coulomb singularities, 
which allows to show, in particular, that the singular term VU which appears in the Liouville 
equation lies in L 1 with respect to the limiting Wigner measure, thereby guaranteeing that the 
weak formulation of the Liouville equation continues to make sense. 

Our methods do not seem to allow to analyse the limit dynamics at eigenvalue crossings, 
since no analogon of our Coulombic non-concentration estimate is available. In fact, the phys- 
ically correct starting point to investigate what happens at crossings would not be QMD, as 
the Born-Oppenheimer approximation underlying QMD also breaks down (see |NW29} IZen32j 
for earliest insights, |Hag94| for a first rigorous account, and e.g. |FG02t ICdV03t ILas04l ILT05] 



for recent results). 

In the remainder of this Introduction we describe our main result precisely. 

Quantum molecular dynamics To simplify matters we assume that all nuclei have equal 
mass. In atomic units (m e = |e| = h — 1), non-relativistic Born-Oppenheimer quantum 
molecular dynamics is given by the time-dependent Schrodinger equation 

'iedt* e (;t)=H e * e (;t) fortGK, 

* e (-,o) = *° (SE) 



e ■ 



with Hamiltonian 



e 2 



H e = --A + U, 



where \P e (-,t) G L 2 (R d ; C) is the wavefunction of the nuclei at time t, d = 3M (M = number 
of nuclei), A is the Laplacian on W 1 , e := (m e /m n ) 1//2 is the (dimensionless) small parameter 
already discussed above (where we have assumed for simplicity that all nuclei have equal mass), 
and U : M. d — > R is the Born-Oppenheimer ground state potential energy surface obtained 
by minimization over electronic states (see e.g. |SQ93j ). The precise definition of U is as 
follows. Let Zi, .., Z M G N and R±, .., Rm £ R 3 denote the charges and positions of the nuclei, 
and let N denote the number of electrons in the system (usually N = Yla=i^ a )- Then for 
x=(R u ...,R M )eR d 

U = E el + V nn , E e£ {x) = inf{^, H ee {x)4j), (2) 



l<a</3<M i=l a=l 1 ' l<i<j<N 

Here the G M 3 denote electronic coordinates and the infimum is taken over the usual subset 
of L 2 ((M 3 x 7j2) n ;C) of normalized, antisymmetric electronic states belonging to the domain 
H 2 ((M 3 x Z 2 ) JV ; C) of H e i(x). Physically, E e £ is the electronic part of the energy, consisting 
of kinetic energy of the electrons, electron-nuclei attraction, and electron repulsion; this part 
depends indirectly on the positions R a of the nuclei since these appear as parameters in the 
electronic Hamiltonian, and can be shown to be bounded and globally Lipschitz (cf. |FriXX] ) 
although it is not elementary to see this. In case N < Yla=i ^a, Zhislin's theorem (see |Fri03] 
for a short proof) says that the infimum in ([2]) is actually attained, the minimum value being 
an isolated eigenvalue of finite multiplicity of H e £. V nn is the direct electrostatic interaction 
energy between the nuclei, and is the origin of the discontinuities of U. 

We remark that quantum molecular dynamics (ISEj) . (pQ), (J2|), which is taken as starting point 
here, itself already constitutes an approximation to full Schrodinger dynamics for electrons 
and nuclei. Its rigorous justification constitutes an interesting problem in its own right; for a 



comprehensive treatment in the case of smooth interactions and absence of electronic eigenvalue 
crossings see |Teu03[ IPST03] . 

The potential (j2j) satisfies the standard Kato-type condition 



U = U b + U s , U b G L°°(R d ), U s (x) = Yl V^Ra-Rp), V a/3 G L 2 (R 3 ) + L°°(R 3 ). (3) 

l<a</3<M 

For such potentials, the operator H e is self-adjoint on L 2 (M a! ) with domain T>(H e ) = H 2 (R d ), cf. 
|Kat51j . By standard results on the unitary group generated by a self-adjoint operator, for any 
initial state ^° G V(H e ) this equation has a unique solution * e G C(R; H 2 (M d ))nC 1 (M; L 2 (R d )), 
the solution operator U e (t) : \P° ^ \l/ e (-,i) being unitary. In particular, 

||tf e (.,t)|| = ||tt°|| forallteR, (4) 
where here and below || • || denotes the L 2 (R d ) norm. 

The Wigner picture Given the state ^> e (-,t) G L 2 (R d ; C) of the system at time t, define the 
associated Wigner function of lengthscale e on R d x R d , 



W £ (x :P: t) = -^— d [ * e (z + f ,t)¥ e (x-f ,t)e-^dy (5) 
(2iT) d J Rd 2 2 

Note that the integrand belongs to L 1 (R^), so W £ is well defined for a.e. x G R d . In fact, the 
integrand is continuous in x with respect to the L 1 (R^)-norm, and hence W e is continuous in 
x. Roughly speaking, W e is a joint position and momentum density of the system. Warning: 
W t is not nonnegative except in the limit e — > 0, but at least its marginals are, 

f 2 

/ W e (x,p,t)dp = ty e (x,i) (position density), (6) 

/ W e (x,p,t)dx = -, -TTT / e~ ip ' x ^ e ty e (x, t) dx (momentum density). 

Jr* {2iie) a / 2 J Rd 

V v ' 

=:(^*«)(PA.t) 

Here and below, denotes the (standard, not scaled) Fourier transform of the function <fi. 
When ty € satisfies (1SED . its Wigner function satisfies 



a,w« = -p ■ v,w, + f„ (we) 

Formally, this follows from a lengthy but elementary calculation which goes back to Wigner, 
see Section H] below. In this section we also introduce a suitable function space setting in which 
the calculation becomes rigorous for general self-adjoint Hamiltonians of type ([1]). We call eq. 
( IWEj) the Wigner equation. Note that it contains no modification of quantum dynamics fISEj) . 
but is just a different mathematical formulation of it. 

Limit dynamics ^From now on we focus on the specific potential energy surface (T5]). In the 
limit e — * 0, the difference quotient in the potential term satisfies 



(due to the Lipschitz continuity of E e £ and the fact that the singular set of V nn is of measure 
zero). To understand what happens with eq. flWEp in the limit, it is useful to split f e into a 
term containing VU(x) ■ y and a term containing the difference quotient of U minus its limit. 
The first term simplifies due to ye~ ipy = iV p e~ ip ' y , giving 



fe = VU{x) ■ V p W e + g £ 



where 

9e = 



(27r) d 



U{X+V V - VU(x) ■ y] + f , t)9 e (x - f , t)e+» dy. 



Formally, passing to the limit in (1WEI) and assuming that g e tends to zero, we obtain the 
Liouville equation 

d t W = -p ■ V X W + W(x) • V P W. (LE) 
This is the transport equation for classical molecular dynamics in R d x R^ with potential U, 

x = p, p=-VU(x). (MD) 

Moreover, assuming the initial data to flSEj) to be normalized, by ([6]) one expects J R2d dW(t) = 1, 
i.e. W(t) should be a probability measure on phase space for all t 6 R. 

The main difficulties in making this rigorous for rough potentials lie in (a) justifying the 
existence of a limiting probability measure on phase space for all times and general initial data 
which are not restricted to 'avoid' the singularities, (b) justifying that g e goes to zero. The 
latter issue arises because when the potential U is not everywhere different iable, the term in 
square brackets does not go to zero for every x, let alone locally uniformly. On the other hand, 
the remaining part of the integrand can concentrate on individual positions x when standard 
semiclassical wave packets such as 

*°Jx) = e- ! 2e iS }- x <t>(2=P), 0<a<l, \\<f>\\ = 1, 



are under consideration, whose Wigner function converges to S( XOtPO ). Thus the only viable 
strategy appears to be to establish that the Wigner function does not charge the singular set 
in the limit. 

Main result Our rigorous result achieves goal (a) in the desired generality, and goal (b) away 
from possible crossings, establishing in particular that the Liouville equation (ILEp remains valid 
across Coulomb singularities. In order to formulate our result we need the following definition. 

Definition 1.1 A sequence of nonnegative Radon measures on R d is called tight if 

lim limsup / dfi e = 0. 



x\>R 



Theorem 1.1 Suppose U : R d — >• R is the Born-Oppenheimer potential energy surface of 
any molecule, or more generally U = £4 + U s with G W 1,oc (R d ), 

U ^= 12 i/^n P C a/3 >0, x = (R u ...,R M )eR d . (7) 



l<a<P<M 



Let {^°} e>0 be a sequence of initial data such that \l/ e G H 2 (R d ), ||^>°|| = I, \\H^° e \\ < c for 
some constant c independent of e, {|^°| 2 } tight Let * e G C(R; H 2 (R d )) n C X (R; L 2 (R d )) be the 
corresponding solutions to the time- dependent Schrddinger equation ( 1SEI) . and let W e be their 
Wigner transforms (JSJ). Then: 

(i) (Compactness) For a subsequence, W t — ^ W in D'(R 2d+1 ). 

(mJ (Existence of a limiting probability measure on phase space) W G C„, ea / C *(R; .M(R 2d )), and 
a probability measure for all t, that is to say W{t) > and J R2 d dW(t) = 1. 

(Hi) (No-concentration estimate at Coulomb singularities) For allt we have VU S G L 1 (diy(t)), 
and W(t)(S x R d ) = 0, where S is the singular set {x = (R%, . . . , Rm) G R 3A/ | R a = R@ for 
some a ^ P with C a /3 ^ 0}. 

(iv) (Limit equation) If Q C R d zs any open sei snca iaai G C 1 (f2), iaen If w a global weak 
solution of the Liouville equation ( 1LEI) on x R d x R, iaai zs to say 



for all(pe C^(Q x R d x R). 

Remarks 1) The regularity requirement [/& G C 1 in (iv) is minimal in order for the weak 
Liouville equation (jSJ) to make sense for general measure valued solutions W. This is due to the 
appearance of Vf/ inside the integral with respect to the measure dW(t). Note however that 
in case of (J2J) this narrowly excludes eigenvalue crossings, as seen from the 2D matrix example 



whose ground state eigenvalue equals —\p\, and is hence Lipschitz but not C . For interesting 
model problems with scalar or vector- valued potentials in which the behaviour of Wigner mea- 
sures past discontinuities of WU can be analysed for suitable classes of initial data see |Ker05] 
and [FG02llEai0lirT05l . 

2) The assumptions on the potential U are far weaker than those needed for uniqueness of the 
Hamiltonian ODE (IMDj) underlying the limit equation. Recall that the standard condition 
guaranteeing uniqueness for ODE's z — f(z) is boundedness of the gradient of the vector field 
/, which in the case of ( IMP I) means boundedness of the second, not the first gradient of U. Our 
assumptions on U are also weaker than those under which uniqueness for weak (L p ) solutions 
to fILEj) is known. The recent nontrivial uniqueness results for transport equations ( |Amb04] ) 
require / G BV, i.e. in case of flMPl) . W G BV (for recent refinements see |BC09j . [AGS08j ). 
Interestingly, however, the latter requirement, while violated by the Coulombic part U s = V nn 
in (j2J), would be met by the model eigenvalue crossing in Remark 1), and expected to be met 
by the electronic part Ub = E e £ in (j2|). We plan to address uniqueness in a separate publication 



3) The higher integrability result in (iii) that VU S G L 1 (dW / (t)), which is essential for making 
sense of the limit equation at Coulomb singularities, requires a quantitative no-concentration 
estimate of form 



for any a ^ f3 with C Q/3 ^ 0. This is because |V^ Q (l/|i? Q - R p \)\ = l/\R a - R p \ 2 > 1/5 2 
in \R a — Rp\ < 5. We do not think that the validity of such an estimate is obvious, the 




8 




|AFFG09j . 




(9) 



naively expected bound only being 0(5) instead of 0(5 2 ) (on grounds of the potential energy 
term J U s \ty e \ 2 , which can be controlled independently of e and t by energy conservation, only 
containing the weaker singularities l/\R a — Rp\). See Sections 4 and 5 for the proof of ©• 

4) In the special case U s — 0, U b G W 2,00 (IR d ), (iv) holds with f2 = M. d , so Theorem 11.11 recovers 
the result of Lions and Paul |LP93[ Theoreme IV. 1.1)]. 

5) In the special case of the potential energy surface (|2J) of the H 2 molecule (M = 2, N = 2, 
Zi — Z 2 — 1), it is known that the ground state eigenvalue of the electronic Hamiltonian is 
nondegenerate. It then follows from a result of Hunziker |Hun86] that Ub is analytic, and in 
particular C 1 , on Q = M 6 \iS, and Theorem [TTTJ justifies classical molecular dynamics globally. 

6) Also, more can be said about the set Q in (iv) in the case of the potential energy surface 
([2]) of a general neutral or positively charged dimer (M = 2, N < Z\ + Z 2 ). By invariance of 
the electronic Hamiltonian H e i and the nuclei-nuclei interaction V nn in ([2]) under simultaneous 
rotation and translation of all particles, we have that U(Ri,R 2 ) = u(\Ri — R2D, that is to 
say the potential is a function of a single parameter, internuclear distance. It then follows 
by combining the result of Hunziker [Hun86] with a classical result of Kato on analyticity 
of eigenvalues of analytic one-parameter families of Hamiltonians [Kat95] that we may take 
fl = M. 6 \(S U C), where S = {R± = R 2 } is the set of Coulomb singularities introduced in part 
(iii) of the theorem, and C = [Jj{\Ri — -R2I = c j} ; the Cj being the (possibly empty) discrete 
subset of 1R + of interatomic distances at which the lowest two eigenvalues of the electronic 
Hamiltonian cross. Note in particular that SUC is a closed set of measure zero; hence Theorem 
11.11 justifies the Liouville equation on an open set of full measure. We expect that when U is 
given by it is always smooth on an open set of full measure. Note however that such a 
result does not follow solely from consideration of non-degenerate eigenvalues as in [Hun86j . 

7) As shown below (Lemma 13. 2p . the weak convergence in (i) also holds in the stronger spaces 
L^ 3 eofc# (R; ^4') and C^ ea / c * i / oc (lR; A'), with A being the Banach space defined in ([To]) . 



2 Wigner-transformed quantum dynamics 

We now make precise in an appropriate function space setting the well known fact (discussed 
informally in the Introduction) that the Wigner transform takes solutions of the Schrodinger 
equation fISEI) to solutions of the Wigner equation flWEI) . Our choice of spaces is convenient 
for our goal to study the limit dynamics for rough potentials. In other contexts other function 
spaces have been considered |Mar89] . 

We begin with the well known formal derivation (assuming that the wavefunction is smooth 
and rapidly decaying for all t). 

Formal derivation Let \l/ e be a solution to (ISEI) . and let W e denote its Wigner transform (jSJ). 
Differentiating the latter with respect to t, we obtain 



d t W e {x,p,t) 



(2k) 



{d t V e (x+ 1, t))V e (x - |, t) + * £ (x + 1, t)d& e (x - f , t) 



-ip-y 



dy. 



By (ISEI) and A± ey / 2 = (4/e 2 )A y this is equivalent to 



(10) 



8 t W e (x,p,t) = f e (x,p,t) 

2i 



+ 



e(2vr) 



{A y V e (x + j,t)) * t (x - f , - *e(x + ^,t)Ay* e {x - ^, - 



QJ 



QJ 



QJ 



-ip-y 



dy (11) 



with f e as in (jWEJ). The Laplacian terms can be simplified via the formula (Aa)b — aAb = 
div(Va ■ b — a ■ V6), an integration by parts, and the formula = ±(e/2)V± ej/ /2, whence the 
integral in ( fTTT) becomes 



div 



ip 



ei 

2i 



-,t))V e (x - —,t) - y e {x + — 



(V B * e (s + f , t)) * e (x - f , t) - V e (x + f ) V„* e (a; - f , f) 



-ip-y 



dy 



{V y V e (x + -f, t)) * c (x -f,t)~ + -f, t)V y V e (x - -f, t) 



-ip-y 



y e (x+^,t)V e (x- e -^,t) 



-ip-y 



dy 



(12) 
dy (13) 
(14) 



(-p-V x W e (x,p,t)). 



Substituting this expression into ffTTj) . we obtain flWEI) . 

Rigorous derivation In the sequel, position coordinates in M. d , d = 3M, are denoted by 
x = (zi, ..,x d ) = (R h .., i? M ) G R d , x< £ M, -Rq e R 3 . 

Lemma 2.1 Tae W%ner transform (TJj o/ any solution ^ e E C(R; H 2 QR d )) n C^QR; L 2 QR d )) o/ 
iae Schrodinger equation flSEj) C/ as m satisfies 



W e E C 1 (R; L°°(R M )), 7^-W e , tt^^™ 

ax,- oxjOx-j 



and solves the Wigner equation kWE]) . 



W e , and f e E C(R; L oc, (M 2d )) /or a// z,j = l,...,d 

(15) 



To obtain an effortless proof, the idea is to express all terms under investigation with the help 
of the following bilinear map which extends the quadratic map \l/ e i— ► W e introduced in (jBI) : 



W*) := 



(2tt) c 



*(x+f)x(x-f)e"^dy. 



Lemma 2.2 Tae map (\l/, x) >— > -^e(^, x) a continuous map from L 2 (R d ) xL 2 (R d ) to L°°(R 2d ). 
In particular, the map \l/ e h- > = -F e (\l/ e , \I/ e ) a continuous map from L 2 (R d ) to L°°(R 2d ). 

Proof of Lemma Q Let x, x' e L 2 (R d ), W = F e (*,x), W = F e (W,x')- Then 

|W(a;,p) - W(x,p)| 

f ((y - + f) x (x-f) + + f)(x-x')(^-f))e- iM dy 

< - + 1)11 llx(* - 1)11 + + 1)11 ll(x - x')(x - f )||) 

= + 11*11 llx-xll) 

for all x and p. Taking \P' = x' — shows W E L°°(R 2ci ), and the estimate above establishes 
the asserted continuity of F e . □ 

Proof of Lemma 12.11 First, we claim that W e E C(R; L°°(R d )). This is immediate from 
^ e E C(R;L 2 (R d )) and LemmaP 



1 

(2n) d 



Next, we investigate the terms -£^W e , aSdx- ^^ anc ^ ^ ne underlying terms 

<9 2 iT/ „„A TT\Ti „™ nflO>. J 2fin>d\\ U„„„„™ iT/ ^ n /"ID) . U2 /TD)d\ \ J + 9 <9 2 



Q^^e and ?7* e are in C(R; L 2 (R d )), because # e G C(R; H 2 (R d )) and the operators 



j dxi ' dxidxj 

and U = H e — yA are continuous maps from H 2 (R d ) to L 2 (R d ). Note now that we have the 
following representations with the help of the bilinear map F e : 



_d_ 

dxi 



dxjdx 



^e = F e (^ e ,^)+F e (* e ,^ e ), 



It now follows from Lemma [2.21 that these terms are in C(R; L°°(R 2d )). 

It remains to show that W e G C 1 (R; L 00 (R 2(i )), with derivative d t W e given by eq. (jWE]). 
First we show continuous differentiability with respect to time. We have 

W e (;t + h)-W e (;t) F £ (V e (;t + h),V c (;t + h))-F e (y e (;t),y e (;t)) 

h h 

Hence by the two convergences M±thhld^l d t W e (-,t) in L 2 (R d ) (from ^ e G C*(R; L 2 (R d ))) 
and # e (-,t + /i) -> tf £ (-,t) in L 2 (R d ) (from ^ e G C(R; L 2 ^))) and the continuity of F e , 

- F e (d t y t (;t),y e (;t)) + F € {* e (;t),dt* e (;t)) in L°°(R M ) as h -> 0. 



Consequently i i— ► W e (-,t) is a differentiate map from R to L°°(R 2a! ), with derivative given by 
eq. (HDD. Continuity in time of the derivative d t W e , i.e. the fact that W e G C X (R; L°°(R M )), 
now follows from ([TOD, ^ e G C(R; L 2 (R d )), d t ^ e G C(R;L 2 (R d )), and - one more time - the 
continuity of F e (Lemma 12. 2p . 

We conclude the proof by showing that W e satisfies ( 1WE|) . The formal derivation of this 
equation from eq. (fTUl) (which we have already established above) has been performed by the 
calculations (jTTj) — (fT4l) . Here, we need only to justify these calculations rigorously. Eq. (TTTT) 
follows immediately from (TTUT) and (1SE[) . Note that all four summands of the integrands on 
the RHS of eq. (ITT]) (cf. also the definition of f e in ( IWEp ). separately belong to L 1 (R^), for 
any x, p and t, because U(x ± f )* e (s ± f , t) and A x ^ £ (x ± f , t) belong to L 2 (R^). Eq. (IT2l) 
follows from the product rule for the Laplacian, and eq. (fT3l) from the fact that the vector 
field inside the square brackets of 012p belongs to the Sobolev space W 1 ' 1 ^) = Wo'lR^) and 
e -ip-y e w 1 ' 00 ^), and that J Rd div(v)(p = - f Rd v ■ V<p for all v G Wj'^R 4 *), <p G W hQoM 



Finally, eq. (I14R follows from an elementary change of variables, concluding the proof. □ 



3 Time-dependent Wigner measures 

We first give a modest technical extension of the construction of Wigner measures in [LP93j . 
Instead of considering a sequence of wavefunctions at a fixed time t, we consider a sequence 
of continuous paths G C(R; L 2 (R rf )) - not required to satisfy any equation - and show that 
under mild conditions these give rise to a continuous path W G C wea fc*(R; A4(R 2d )) of Wigner 



measures. See Lemma 13 .21 We then combine the lemma with careful a priori estimates on the 
singular contributions to f e in eq. (IWEI) and a lemma on propagation of tightness under (jSEft 
to prove Theorem 11.11 (i) and (ii). An interesting feature of these proofs is that, in contrast 
to the existing literature, they are extracted directly from the Schrodinger dynamics, without 
relying on a representation of the limit measure as push-forward of its initial data under (IMDj) . 
which is not available here. As regards our proof in (ii) that f R2d dW(t) = 1 for all t, our 
argument via Schrodinger dynamics is inspired by Corollary 1 in [FL03] (see also Proposition 
4 in [Las04| ). 

We begin by recalling the notion of weak convergence of Wigner transforms at a fixed 
time t to Wigner measures introduced in |LP93j . For an alternative construction of Wigner 
measures see |Ger91aj . Closely related constructions are the microlocal defect measures of 
Gerard [Ger91b] , motivated by questions in microlocal analysis, and the if- measures of Tartar 
|Tar90j . motivated by questions in homogenization theory. 

Let A denote the following Banach space 

A :={(/) e C (IR 2ci ) | ||0|U := / sup|(^)(x,y)|dy<oo}. (16) 

JR d x£R d 

Here Co(M 2a! ) is the usual space of continuous functions on M. 2d tending to zero at infinity, and 
JF p is the partial Fourier transform (jF p 0)(x, y) = J Rd e~ zp ' y 4>(x,p) dp. Since A is a dense subset 
of C (M 2a! ), its dual A' contains Co(M 2d ) = M(M. 2d ), the space of not necessarily nonnegative 
Radon measures on M. 2d of finite mass. In particular, the delta function o~( XQtPo ) centered at a 
single point (xo,Po) in classical phase space belongs to A', and weak* convergence in A' allows 
convergence of smeared-out Wigner functions coming from a quantum state to a delta function 
on phase space (i.e., a "classical" state). 

The following basic facts were established by Lions and Paul. 

Lemma 3.1 (Wigner measures) lLP93f 

(i) (compactness) Let {^ e } be any sequence in L 2 (R d ) such that 

\\<!> e \\ 2 <C (17) 

for some constant C independent of e. Then the sequence of corresponding Wigner transforms 
{W e } contains a subsequence {W t >} converging weak* in A' to some W G A' . 

(ii) (positivity) Any such limit W G A' satisfies W G Ai(R 2d ), W > 0. In other words, W is 
a nonnegative Radon measure of finite mass. 

(Hi) (upper bound) Let {^ e '} be a further subsequence such that \^ e '\ 2 /i in .M(IR d ). Then 

[ W(;d P )<^. 

In particular, any limit W as in (i) satisfies j R2 d dW < C , with C as in (T7\ ). 

(iv) (preservation of mass) If H^H 2 = C for all e, and the sequences of position and momentum 
densities are both tight, that is to say 

limsup / |\l/ e (x)| 2 dx — > 0, limsup / 

e^O J\x\>R e^O J\p\>R 



(2ire) d /- 



V 



dp -> (R -»• oo) 



then / dW = C. In particular, in the case C — 1, W is a probability measure. 
Jm: m 



For future purposes we note that (i) is immediate from the Banach-Alaoglu theorem and the 
elementary estimate 

W^d(x,p)\<j^\\^ e \\ 2 U\\ A for0eA, (18) 



which implies 

||W e |U'= sup < * ||* e || 2 , (19) 

teAm \m\A (2<r 

that is to say {W e } is a bounded sequence in A'. The proofs of (ii) and (iii) are less elementary, 
and require use of the Husimi transform; see [LP93j . 

An analogue yielding time-continuous paths of Wigner measures is the following. 

Lemma 3.2 (Time- dependent Wigner measures) 

(i) (compactness) Let {^ e } be a sequence in C(M;L 2 (IR d )) such that 

sup||v]/ e (.,t)|| 2 <C (20) 

teR 

for some constant independent ofe, and let {W e } be the sequence of associated Wigner functions. 
Then for a subsequence, W e — W weak* in L°°(M; A'). 

(ii) Suppose in addition that for any test function G Cg°(M 2d ) the functions 



feAtf : = W e (x,p,t)(j)(x,p)d(x,p) 
JR 2d 

are differentiate and satisfy 

sMitkM <c* (2i) 

teR 

for some constant C$ independent ofe. Then, W G C wea ^(M.] Ai(M. 2d )), W(t) > for allt, and 
W e (-,t) — Wit) in A' for all t. Moreover, the latter convergence is uniform on compact time 
intervals, i.e., for any test function G A and any compact /cR, fe,<f>{^) converges uniformly 
with respect to t E I to J R2d <f> dW(t). 

Here (i) is a straightforward adaptation of the time- independent theory in [LP93j . The key 
point is the assertion in (ii) that the limit measure has slightly higher regularity in time than 
naively expected (continuous instead of L°°). This allows, in particular, to make sense of initial 
values. 



Proof The first part is an easy consequence of Lemma I3TT1 cf. ( 1191) . and assumption (1201 . which 
imply that {W € } is bounded in L°°(1R; *4/). Since the latter is the dual of the separable space 
L 1 (M;^4), the assertion follows from the Banach-Alaoglu theorem. 

The second part is less trivial, and requires various approximation arguments. First, we 
test the weak* convergence of {W e } from (i) against tensor products 4>(x,p)x(t) with <fi G A, 
X G L 1 (IR). This gives 



f e At)x(t)dt= W e cj)®xd(x,p,t)^ (/>® X dW(t)dt= I U(t) X (t)dt, 

jR?d+l J R J R 2d 

where f^t) := J R2d (f> dW (t) . Consequently f e><t> A in L c 



2d ). Then by assumption (I2ip and the compact embedding W 1,00 ([— T, T]) 
converges uniformly to on any compact interval [— T, T] . In particular 



Now let 
C([-T,T]), '/, 

uniform limit of continuous functions, continuous, and f e ,(/>(t) 
t G R, that is to say 



W 6 (-,t) W(t) in £>'(M M ) pointwise for all t G 



/</>(*) pointwise for all 



(22) 



Now fix t. For a further subsequence which may depend on t, by Lemma [3.11 W e /(-, t) — ^ 
in .4' C £>'. Together with ([22]) this yields W(t) = W(t), as well as the convergence W e (-,t) 
W(t) in A' for the whole sequence. By Lemma 1371] all remaining statements about W(t) follow, 
except its asserted continuity in t. 
It remains to show the latter, i.e. 



<dW(t + h) 



<dW(t) (h-*0) forall^GC 



>2d\ 



Given G Co(M 2d ) and 5 > 0, by the density of C~(M M ) in C (l 
such that Wdfj — < $, where II • denotes the norm of L°°( 



l 2d ) there exists (ps 
¥L n ). Consequently 



GC-( 



(23) 



!>2d\ 



< 



< 



<j>dW(t + h) - 
4> s dW(t + h) 



<j>dW(t) 



</> S dW(t) 



(4> - 4> s ) dW(t + h) 



4>s) dW(t) 



[ <f> s dW{t + h)- I <f> s dW{t) +6( f dW(t + h)+ [ dW(t)). 

JU 2d JR 2d VR2<* JR 2d ' 



As h — > 0, the first term vanishes by the continuity of J" R2d 4>&dW{t) = f<fi s (t) in t (which was 
already established above). The second term stays bounded by 2C5 by Lemma I37T1 (iii), giving 



lim sup 



>dW(t + h) 



•dW(t) 



< 2C5. 



Since 5 was arbitrary, the continuity assertion ([231 follows, and the proof of Lemma 13.21 is 
complete. □ 

Remark Functional analytically, Lemma [3721 says that under the assumptions f[20l) . fl2TT) . {W e } 
is relatively compact in C wea k*,eoc{^-] A'). This may be viewed as a weak-convergence variant 
of the well known compactness lemma of J. L. Lions |Lio69[ Chap. 1, Theoreme 5.1], in which 
the condition of boundedness of time derivatives in some Banach space has been replaced by 
condition (|2~TT) which is related to a weak topology. 



We close this section by applying the above lemma to prove the first two statements of Theorem 

o 

Proof of Theorem Q7J] (i) This follows from ||^ e (t)|| = H^H = 1, Lemma [372] (i), and the 
fact that weak* convergence in L°°(IR; V 4 / ) implies convergence in V(R 2d+1 ). □ 

Proof of Theorem Q7J] (ii) First we prove that W G C weak *(R; M(M? d )), W(t) > 0. By the 
Lemma 13.21 all we need to show is that for any function G C[f(IR M ) on phase space, the 
expected value 

ferf(t) = W e (x,p,t)(/)(x,p)d(x,p) 
jR2 d 



is differentiable in t and satisfies hypothesis fl2T|) of Lemma 13721 (ii), i.e. f' e At) stays bounded 
independently of e and t. Differentiability in t holds even for £ L 1 (M 2c( ), since W e £ 
C 1 (M;L°°(M 2d )) by Lemma EU 

To deduce ( 12TT) we start by exploiting the Wigner equation (IWEjl . This yields 



\f'^(t)\ < I W e (x,p,t)p- V x <j)(x,p)d(x,p) 

U b (x + f ) - U b (x - f 



+ 
+ 



(2?r) d 
1 



ey ev 
■V e {x + y , t)tf e (s - -f, t){T p <f>){x, y) d(x, y) 



(27T) 



f/ s (x+f)-f/,(x-f) 



# e (x + ^, t)tt 6 (z - ^, t){F p <f>){x, y) d{x, y) 



The first term on the right hand side is bounded by -^^\\^/ t (-,t) || 2 ||p-V x ^||^i by (JTSJ), and hence 
bounded independently of e and t, by (jlj) and = 1. 

Thanks to the elementary inequality | — -\ < || V£4||oo|?/|, the second term is 

bounded independently of e and t by 



(2tt) c 



IW, 



6 



SU P l(^>0)(^,2/)l dy. 



Finally, in order to estimate the third term we observe that for R, Q £ W 1 



R 



R 



I 



< 



\Q\ 



\ R+ ¥1 \ R 



eQ\' 



Hence, with x = (Ri, . . . , Rm), y = (Qi, • • • , Qm) and setting R = R a — R/3, Q = Q a — Q/3 £ M 3 , 
the third term is bounded by 

TTW E IP 1 B | ^(-»*)H 2 / IQa-Q^SUpK^^^ldy 



l<a<P<M 
(27r) d 771 



1 2/ 1 sup |(.7y/>) (x, 2/) | dy, 



where m = min{C aj g | C a p ^ 0}. But the right hand side stays bounded independently of e 



and t thanks to Lemma 15.11 establishing hypothesis (1211) and thus completing the proof that 
W £ C weak *(R;M(M. 2d )), W{t) > 0. 

It remains to show that J R2d dW(t) = 1 for every t. In Lemma [3.21 (ii) we proved that W £ (-,t) 
converges weak* in A' to W(i) for every t. Hence Lemma 13.11 (iv) is applicable and it suf- 
fices to verify that both the sequence of position densities |\l/ e (-,t)| 2 and momentum densities 
| ^ 27re ^d/2 (-T^eXf ; t)\ 2 are fight. As regards the momentum densities, this follows from uniform 
boundedness of kinetic energy, (153!) . since 



p\>R 



(27re) d / ; 



1 1 

dp < — -- 

y ~ R 2 2 



~ [ |eW e (x,t)| 2 dx 
R 2 J Rd 



(W e )(Z,t) 2 \p\ 2 dp 



Finally, tightness of the position densities follows (under much weaker hypotheses on potential 
and initial data) from the Lemma [3.31 below, completing the proof of Theorem 11.11 (ii). □ 



Lemma 3.3 (Propagation of tightness) Let U be as in ([3]) , and suppose in addition U bounded 
from below. Let {^ t } be a sequence of solutions to the time-dependent Schrddinger equation 
(Elf, whose initial data satisfy \&° £ H 2 (R d ), \\^° e \\ = 1 (normalization), and H e W°) < 
C (bounded energy). If the sequence of initial position densities {|^°| 2 } is tight, then so is 
{|^ e (-,t)| 2 } ; for allte R. 

Proof of Lemma [331 Let x e Cg°(R d ), < x < 1, X = 1 on \x\ > 1, x = on |z| < 1/2, 
and for i? > set := x( x l R)- Then |Vx_r| < C/i2, |A%_r| < C/i? 2 for some constant C 

independent of R. Letting (•, •) be the L 2 (R d ) inner product and abbreviating (A)^ = (0, A<f>), 
we obtain from flSEj) 

d i 

Since [xii, fl" e ] = y A Xr + e 2 V^i? • V, 

4<X*W,t) < ^ / |Axfl(x)| |^(x,t)| 2 dx + / |eV^(x,t)||Vxi?(x)||^(x,t)|dx 
< | IIAxfllU ||* 6 (-,<)ll + IIVxhIIoo ||eW e (-,t)|| ||* e (-,t)||. 

^From the boundedness from below of U and the conservation in time of the energy (^ r e (-, t), H e ty e (-, t)), 
we obtain ||eV 1 S r e (-, t)\\ < const, for some constant independent of e and t. Using the bounds 
on IIVxrIIoo and HAxrHoo, and assuming without loss of generality e < R, it follows that 
■§t(XR)^e(-,t) < c °R St ' f° r some constant independent of e and t. Consequently (considering 
without loss of generality t > 0) 

/ |^(x,t)| 2 dx< (Xfl)*«(.,t) < / —^~ds< |^( x )| 2 dx + — — t^O 

J\x\>R Js=0 K J\x\>§ K 

as i? — > oo, by the tightness of the sequence of initial position densitites. □ 



4 Justification of the Liouville equation 

We now show that time- dependent Wigner measures arising as a limit of semiclassically scaled 
solutions to the Schrddinger equation are weak solutions of the Liouville equation (1LED . 

In line with the discussion in the Introduction, the main work will go into analyzing the 
behaviour near the Coulomb singularities. The first order of business will be to verify that the 
limit equation even makes sense, for which we need that the gradient of the potential lies in L 1 
with respect to the Wigner measure, for all t. 

We note the particularly impeding feature that the singularities of the Coulomb forces 
appearing in the Liouville equation have the magnitude l/\R a — Rpl 2 , and are hence much 
worse than the singularities of the Coulomb potentials l/\R a — Rp\ appearing in the original 
Schrddinger equation. 

Proof of Theorem 11.11 (iii) We want to show that 

^)(5xl d ) = 0, (24) 
/ |W S | dW(t) < oo for all t £ R. (25) 



To this end, we fix t G K. and choose a subsequence such that |W e (-, t)\ 2 — ^ in ^(M '). 
According to Lemma 13.11 (iii) we have 

/ W(-,dp,t) </i. (26) 
Jp<=m. d 

First we show W(t)(S x R d ) =0. Bv (I2B|) it suffices to show that /x(5) = 0. Let G C(M d ) 
with < < 1, = 1 for dist(x, S) < 5, and = for dist(x, S) > 25. We use the elementary 
property of weak* convergence in Ai(M. d ) that if /j, e > 0, fj, e — fi in Ai(R. d ), and / is a bounded 
nonnegative continuous function, then J Rd f dfi < liminf J Rd f dfi t . (This follows immediately 

by approximating / by compactly supported functions fj with < fj < f and dominated 
convergence.) Consequently 

/ < / 0d/i< liminf / |\P e (-, t)\ 2 dx = liminf / |\P £ (-, t)| 2 dx. 

Jdist(x,S)<S JR d JR d Jdist(x,5)<2<5 

The idea now is to use the fact that U s > on {x G M d | dist(x, S) < 25} for some constant 
C, and appeal to Lemma IBTTl below. This yields 

/ dfi < te) 2 liminf / U 2 \^ e (-,t)\ 2 dx < (csYc. (27) 

Jdist(x,S)<5 ^ ' JR d ^ ' 

Since the RHS tends to zero as 5^0, fi(S) = and W(t)(S x R d ) = 0, completing the proof 
of (JMD. 

To show fl25l) we will use the monotone convergence theorem. To this end we set for 5 > 



f 5 :=mm[f,^Y f := |W S |. 



Then is a bounded continuous function on M. d . By the fact that \VU S \ < C$U 2 (thanks to 
the identity V R a ^ = ^ J^p ) an d Lemma [57TI below. 

/ / 5 d/i < liminf / / a |* e (-, t)| 2 dx < C liminf / U 2 \^ t (-,t)\ 2 dx < C C, (28) 

where Co, C are constants independent of e and t. 

Consider now the limit 5^0. In this limit, fs —+ f monotonically on Since /J*(S) = 0, 

it follows that fg — > f /i-almost everywhere. Hence the monotone convergence theorem yields 
/ G L^dyu) and 

fdfi = hm [ f s dfi. (29) 



Consequently, by the p-independence of /, fl26|) . fl29l) . (|28|) . 

/ / dW(t) < [ fdfi<C C<oo, (30) 

establishing (1251) and completing the proof of (iii). □ 

Proof of Theorem 11.11 (iv) Let f2 C W 1 be an open set such that Ub G C 1 (f2). We need to 
show that eq. ([8]) holds for every G Cq° (f2xIR d xlf). Starting point is the fact that by Lemma 



12. H W € satisfies the Wigner equation (jWEft . Multiplying flWEp by 4> and integrating by parts, 
we obtain the weak form 



/ (w e (dt + p-V x )<l> + f e <f)d(x,p,t)=0. (31) 

Passage to the limit e — * in eq. ([3"Tj) is done in three steps, carried out in the order of 
increasing difficulty: 1. Analysis of the local terms, 2. Analysis of the nonlocal term J f e <p for 
test functions vanishing in a neighbourhood of the Coulomb singularities, 3. Analysis of the 
nonlocal term in a neighbourhood of the Coulomb singularities. 

Step 1: Analysis of the local terms in (I3TI) 

As e — > 0, W e W in L°°(1R; A'), and hence in particular in T>'(R 2d+1 ). Thus, the first term in 
(l3"Ij) satisfies 



W e (d t + p-V x )<f>d(x,p,t)^ / (d t +p-V x )(f>dW(t)dt (e->0). (32) 




Step 2: Analysis of the nonlocal term J f t <p for test functions vanishing in a neigh- 
bourhood of the Coulomb singularities 

Let S C R d be the set of Coulomb singularities of the potential U s (see Theorem ll.il (in)). Our 
goal in this step is to prove that 



/ f e <f>d(x,p,t) -> - / / VU -V p( f)dW{t)dt (e->0) (33) 



for all test functions on f2 which vanish in a neighbourhood of S, i.e., G C§°((f2\iS) x 

R^xl). This together with (J3T1) . ( J32l) completes the proof of (iv) for the above test functions. 
We begin by rewriting the left hand side of (1331 . Substituting the definition of f e and 



carrying out the integration over p gives (abbreviating = ty e (x + y,t), ^ e = ^e{x — ^f,t), 
T p (j) = (J" p 0)(x,y,t)) 

r j r jj( x _l m) - JJ(x - 

/ f e( pd(x,p,t) = -—- d -± 1± ^ ^ e * e (f p <j>)d(x,y,t). (34) 

The idea now is to split the ^-integration into two regions in such a way that either \ey\ << 1, 
in which case the difference quotient of U is well approximated by the derivative VI7(x) ■ y, 
or \y\ » 1, in which case T p 4> is very small, due to the rapid decay of the Schwartz function 
as \y\ — > oo. To implement this idea, we introduce a cut off radius e~ a with some fixed 
a G (0, 1) and choose the regions of integration as \y\ < e~ a and \y\ > e~ a . In particular, 
we denote A e := {(y,t) G R d+1 \ \y\ < e~ a } and V e := {(y,t) G R d+1 \ \y\ > e~ a }. It is also 
convenient to subtract off, and then add again, the RHS of (1341) with the difference quotient of 
U replaced by the derivative VU(x) ■ y. This yields the following natural splitting of the RHS 



of ([34]) into a sum of five terms 



/ 



1 2 



T 3 
T 4 



f e <j>d(x,p,t) = T x + T+ + T 2 ~ + T 3 + T 4 
-f/(x + f)-f/(x-f) 



(35) 



(271")' . 

% 



1 



- W(x) • yj J-, • |y| d(x, i), 



=F 



(2vr) c 



U(x±f) — 



* e \P e (^)d(a;,2/,t), 



(27r) d j M d x 



(2tt) c 



W(x) -2/* e W 6 (^) d(x,y,t), 
W(x) ■ y yWe{Fp<t>) d(x,y,t). 



Here, in Ti the factor rr • lyl = 1 has been inserted for future use. 

' 1 \y\ m 

We begin by analysing T±. Let 5, R > such that the set := {x G | |x| < 

R, dist(x, S U 9f2) > 5} contains all x with (x,p,t) G supp 0. Then provided e is sufficiently 
small 

\x±^-\<2R, dist(x± ^,5Uffi) > - for all x G \y\ < e" Q , 

that is to say x ± y G S^2_r,<5/2- Hence by the continuous differentiability of U in ^2R,s/2, 



U{x + f ) - - 



-) y 

► VU(x) ■ — as e — > uniformly for x G Oflj, lyl < e 

1 2/ ] 



Consequently, applying the Cauchy-Schwarz inequality with respect to the integration over x, 
T\ can be estimated by 



\T X \ < j—rrg SUp 

l^rj xen R , s ,\y\<e- 



\y\ 



£\y\ 



sup ||# e (-,t)|| 2 \\yF P <f>\\*^ o 



as e — > 0, where here and below, for any function x G 5(M 2d+1 ) we denote 



sup |x(z,y,t)|d(y,t). 



(36) 



(37) 



The terms can be estimated in an analogous manner, again applying the Cauchy-Schwarz 
inequality with respect to the integration over x: 



\T?\ < 



(27T) 



sup \\U^ e 



sup||tt £ (-,t)||- / sup\(F p <l>)(x,y,t)\d(y,t). (3i 



By Lemma [57T1 and the boundedness of Ub, the norm ||t/" 1 S r e (-, t)\\ stays bounded independently 
of t and e. On the other hand, since g 5(IR 2d+1 ), 

sup (l+\y\) m \(F P 4>)(x,y,t)\ =: c m < oo 

(x,y,t)m 2d+1 

for any m — 0, 1, 2, . . .. Consequently, for all m > d + 1 and all |y| > e~ a 



\(F p <f>)(x,y,t)\ < 



< 



c a(m—d—l) 



(l+|y|) m ~ (l+\y\) d+1 ( e -<*)m-<i-i (l+|y|) 



d+l 



Hence, choosing T so large that |£| < T for all (x,p, t) G supp <j), the last factor in (138]) satisfies 



sup I (F p cj>) (x, y, t) | d(y , t) < e"^-^ 1 ) / " dy . (39) 

v e xm d Jm d \ L +\y\) 



Then, if m is chosen so large that the exponent a(m — d — 1) > 1, the positive power of e in 
(|3~9~1) 'beats' the singular factor i in (138]) . Thus 

if -> as e -»• (40) 

(in fact, faster than any power of e, but this is not needed in the sequel). 

The proof that T 3 — > is analogous but easier, due to the absence of a singular prefactor - 
and the fact that Vt/(x), unlike U(x ± ey/2), is bounded on supp^. We simply estimate 

|T 3 |<— ^ sup |W(x)| [ \y\sup\{F p <j>){x,y,t)\d(y,t)^0ase^0 (41) 

by P2J with m> d + 2. 

Finally, consider the last term, T 4 . Using iyT v $ = ^ P (V P 0) and interchanging the integra- 
tions over p and y, 

T 4 = --^ / VU(x)-\[ V p( f>(x,p,t)e- ip - v d P U e (x + ^,t)y e (x-^-,t)d(x,y,t) 
(27r) fl J R 2 d+ i LJjgti J 2 2 

W (x) • V p 0(x, p, t) W e (z, p, t) d(x, p, f). 



Since Vf/ ■ V p G L 1 (M; .4.) due to the continuity of V£7 on supp0, and since W e — ^ W in the 
dual L°°(1R; ^4'), 

T 4 -»■ - / / W(x) ■ V p 0(x, p, t) dW(f) dt as e -> 0. (42) 

JR JM. 2d 

Combining ([36]), gDD, (3D, (H3) yields (133]) for G C~((fi \ 5) xR-'xR). 

Step 3: Analysis of the nonlocal term J _f e in a neighbourhood of the Coulomb 
singularities 

We prove here that eq. (133]) continues to hold for arbitrary (p G C^°(fi x M. d x M.) not required 
to vanish in a neighbourhood of the Coulomb singularities, i.e. that the Liouville equation 
continues to hold across Coulomb singularities. 

This is quite remarkable, since the available a priori bound f Rd U 2 |^ r e (-, t)\ 2 < const, only 
rules out concentration of the measure |^ r e (-,t)| 2 on Coulomb singularities (as was shown in 
(]27j)). but not concentration of the blown-up measure U 2 \^ e (-, t)\ 2 ~ |VZ7 S | t)\ 2 which 
asymptotically appears in f f t (f> (see the leading term T 4 in (135])). This suggests the possibil- 
ity that an additional contribution of form f R fs x m. d $ dv(t) dt, with u{t) a singular measure 
supported on S x R d , could appear in the limit equation (jSJ). 

This possibility will be ruled out by careful use of the evolution equation (jWEI) satisfied by 
W t . Roughly, our analysis below will lead to the insight that the asymptotic amount of mass 
of f e in a (^-neighbourhood of the singular set S is at most of order 5, not order one. 

We will need rather precisely chosen cutoff functions. Given 5 > 0, we let 

\R n — R 



T] S {x) := 

l<a<p<M 



II »/( a ) with * = 



where 77 G C^°(R) with < 77 < 1, 77 — 1 on \z\ < 1/2, rj = on \z\ > 1. Then for some 
constants Ci, C 2 independent of 5 



Q 

iVrjsix) I < — for all x, 775 (x) = for distfcc, <S) > C2S. 




(43) 



We now write = (1 — 775)0 + 77,50, and consider both contributions to J R2d+1 f e (j) d(x,p,t) 
separately. Since (1 - rj s )(j) belongs to Cg°((fi \5)xR rf xl),we have by Step 2 (cf. (133"]) ) 

f e (l- Vs )<f>d(x : p } t)^- [ [ VU-(l-7i s )V p <t>dW(t)dt ase^O. (44) 

We now claim that the remaining terms are small when 5 is small, that is to say 



lim sup 

e-*0 



f e T] S <f)d{x,p,t) 



as 5 — > 0, 




VU ■ rj S V p ct)dW{t) dt -> as 5 -> 0. 



(45) 
(46) 



Clearly, (1341) and (1431). (HBj) imply eq. ([33]) for the arbitrary test function G Cg°(fi x R d x R). 

Proving (146]) is not difficult, but for convenience of the reader we include a proof. Consider 
first a fixed t. By the facts that < 775 < 1 and supp7/5 C {x G R d | dist(x, S) < C 2 S}, 



VU ■ m V p ct>dW{t) < ||V p 0|U / |W| dW(t) (47) 

J (supp r^nsupp <j>) xR d 



\VU\dW(t) = as 5 -> 0, 

'(cSnsupp </>)xR d 

since dW(f)(S x R d ) = and V£7 G L 1 (supp0; dW(f)). Moreover the LHS of (JSJ) stays 
bounded independently of t by (1301) . Hence the integrand with respect to t in (|46|) tends to 
zero boundedly a.e. as 5 — -> 0, and so by dominated convergence we infer (|46|) . 

It remains to establish (j4"5l) . This is the difficult part of Step 3, due to the fact discussed 
above that the a priori bound of Lemma 15.11 does not rule out the possibility of concentration 
of mass of | Vf/||^ r e (-, t)\ 2 on the set S of Coulomb singularities. Using first (jWEp and then the 
definition of W f we have 



f e r] S (j)d(x,p J t) 
W, 



W e (d t +p-V x )(r) S 4>)d(x,p,t) 

1 



(27r) d J^d xAC 



Vs(d t +p-V x )(f) + Vr]s-p(f) d(x,p,t) 
1 



(2k) 



u eiS d(x,y,t) 



u £)S d(x,y,t) 



(2k) 



u etS d(x,y,t), 



=:Ql 



[ xV c 



where we have split the domain of integration as in Step 2, and where 



u e ,s(x, y, t) = r)s(x) T p ((d t + p ■ V x )(f>)(x, y, t)+Vr] S (x)-J r p (p(p)(x, y, t) ty t (x+— , t)V e (x - — , t). 

s „ ' v — — ' J z z 



The second term, Q 2 , is exactly of the same form as the term T 3 in Step 2, with VC/(a?) 
replaced by r]$(x) respectively Vrjs(x), and the test function yT p § replaced by x respectively 
£. Consequently, estimating as for (14T1) . 



IQ 2 | < 



(27T) 



sup |x(z,y,*)|d(y,*) + l|V7fe||oo / sup d(y,t) -»0(e -►()), 



V t x 



(48) 

by (159~j) with x respectively £ in place of ^,0, and m > d+ 1. The first term, Qi, will be dealt 
with by an argument similar to the no-concentration estimate 02 7|) in the proof of (iii), except 
now a non-local version is needed since the term under investigation is not local in W e . For e 
sufficiently small, and (x,y) belonging to the domain of integration, that is to say x G suppr^, 
\y\ < e~ a , we have (writing x = (R u .., R M ), V = (Qi, Qm)) 

\{R a ± e ^)-(R p ± e -9^)\<25 forl<a</3<M 



and consequently U s (x ± y) > ^ f° r x e SUPP 7 ?^ M < e a ; an d some constant C. Hence 



IQil < 



< 



(C5f 
(27r) d 

(C^ 
(27r) d 



! xA e 



Vs\x\ + |V^| |e|] |E/ s (x+^)tt £ (x+^,t)| |f/ s (x-^)^(x-^,t)|d(x,y,t) 



MUlxll* + llV^Hoolieil, sup \\U s V £ (;t)\\ : 



with || ■ ||* given by fl37j) . Since ||f/ s ^ r e (-, £)|| stays bounded independently of t and e by Lemma 
15. H Halloo = 1, and HV775H00 < ^ by fT4"3]) . it follows that 



\Qi\ < C*[5 2 + 5} 



(49) 



for some constant C* independent of e and 5. Note how the positive power of 5 gained by 
inserting the multiplier U s has 'beaten' the negative power of 5 coming from the gradient of 
the cutoff function rjg- Combining (I48p . 049]) gives 045 p . This completes the proof of Theorem 
11.11 (iv) for general test functions. □ 



5 An a priori estimate for the Schrodinger equation with 
repulsive Coulomb interactions 

We prove now the a priori estimate used in the derivation of the Liouville equation that the 
potential term U s ^> e (-,t) in the semiclassically scaled Schrodinger equation flSEj) stays bounded 
in L 2 (R d ) independently of e and t. 

For potentials with Coulomb singularities, such as (J2J), such an estimate says in particular 
that in the limit e — > 0, the wavefunction cannot concentrate mass at the singularities. 

On physical grounds, one would expect this to be true only for repulsive interactions (as 
present here), but not for attractive interactions. The challenge then is to translate this phys- 
ical intuition into a mathematical argument fine enough to detect sign information. This is 
achieved by the positive commutator argument below, which exploits not just the repulsivity 
(i.e., positivity) of U s , but also its special Coulombic nature. 



Lemma 5.1 Let U = U b + U s be as in Theorem\TJ\ Let {\1>°} C V{H e ) satisfy ||^°|| = 1 and 
ll-^e^ell — c for some constant c independent of e. Then the solution ty e (-,t) to ( ISEj) satisfies 



sup \\U a V e 



< C 



(50) 



for some constant C independent of e. 

Proof By standard results on the unitary propagator e~ ltH associated to a self-adjoint operator 
H, if e V(H), then so is ip(-,t) = e"^ for all t e R, and \\^(-,t)\\, (tp(-,t), Hip (■,{)), 
\\Hip(-, t) || are time-independent. (Formally, the time-independence follows from the Heisenberg 
evolution equation for expected values, ■^(ij)(-,t),Aip(-,i)) = (ip(-,t),l[A, H]ip(-,t)), by taking 
A = I, H, H 2 .) Applied to our case this yields, besides (@J, 



sup(- / |eW £ (x,t)| 2 dx + / (U b (x) + U s (x))\^ e {x,t)\ 2 dx) < const. 
teR ^2 J R d J R d J 



sup 



< const.. 



(51) 



(52) 



the constants being independent of e. By the boundedness of [/&, and the nonnegativity of 

U s , 

sup- / |eV^ e (a;, t)\ 2 dx < const., sup / U s (x)\^f e (x, t)\ 2 dx < const., (53) 

tm 2 J R d ten Js.d 

all constants being independent of e. Now we expand the left hand side of ( 1521) . and rewrite 
the latter in the form 

sup( (~A + U b )y e (-,t) 2 + 2Re(-^Atf £ (-,t),£/ fl * £ (-,t)\ 

+ 2Re(f/^ £ (^),^ £ (-,t)) + ||^ £ (-,t)|| 2 ) < const. (54) 

Using the positivity of U s and ( l53l) . the third term satisfies 

sup|2Re<W e (-,t),W 6 (-,t))| < 2||C/ b || 00 sup|(^ e (-,t), ^ £ (-,t))| < const. (55) 

teK ten 

The key point now is the following claim: 

Re(-A^, U s ip) > for ^ e H 2 (M d ). (56) 
Postponing its proof, substitution of (jo"5]) . (1561) (with ^ = \l/ £ (-,t)) into (1541) yields 



sup 

teM. 



A + U b )ty e (-,t) < const., sup \\U s ^ t (-,t)\\ 2 < const., 



establishing the assertion. 

It remains to prove (|56|) . (This depends on the Coulombic nature of U s as well as the fact 
that it is positive, i.e., repulsive.) By a standard approximation argument, using the density 
in H 2 (M ci ) and the fact that, by Hardy's inequality, ip i— > U s ip is a continuous map 



of Cg° 



from H 1 (R d ) to L 2 (IR d ), it suffices to prove (J5| for if) e C 



In this case, compute 



Re(-A^, U s i[)) = Re / ■ V(tf» 



|V^| 2 ^ + Re / (WOV- VC/, 



The first term is > and the second term equals 

If ((V^V + W).W a = I [ V\^\ 2 -VU s = \f |V| 2 (-Af/ s ). 

* JR d * JR d * JR d 

Now, considering e.g. the term in U s , cf. (j7|), write 

\R± — -R2I 

— A(Rl,...,R M ) = -A fli+Hg - A Rl -R 2 - A : /,. ;i /,.,.: 

and hence 

-A^,,,,,^) 1 = -2A fil _ fl2 - 1 = 87r5(i? 4 - i? 2 ) 
l-Ki - n 2 \ \tti - n 2 \ 

(recall -AA = 4tt(5 in IR 3 ). Thus / Rd |^| 2 (-Af/ s ) > 0, which completes the proof of §§]). □ 
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